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Abstract 

An anisotropically forced passive vector model is analyzed at scales much smaller and larger than 
the forcing scale by solving exactly the equation for the pair correlation function. The model covers 
the cases of magnetohydrodynamic turbulence, the linear pressure model and the linearized Navier- 
Stokes equations by choice of a simple parameter. We determine whether or not the anisotropic 
injection mechanism induces dominance of the anisotropic effects at the asymptotic scaling regimes. 
We also show that under very broad conditions, both scaling regimes exhibit anomalous scaling 
due to the existence of nontrivial zero modes. 
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I. INTRODUCTION 



One of the most important problems of turbulence is the observed deviation of Kol- 
mogorov scaling in the structure functions of the randomly stirred Navier-Stokes equations 
in the inertial range of scales [l^]. Contemporary research in turbulence has recently 
provided an explanation for this phenomenon in the context of passive advection models 
(see e.g. [l^ for an introduction and further references). In the case of the passive scalar 
model describing the behavior of a dye concentration in a turbulent fluid, such a violation 
of canonical scaling behavior (henceforth referred to as anomalous scaling) has recently 
been traced to the existence of a type of statistical integrals of motion known as zero 



modes 
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The result can be obtained under some simplifying assumptions about the 
velocity field, namely assuming the velocity statistics to be gaussian and white noise in 
time, which results in a solvable hierarchy of Hopf equations for the correlation functions. 



Such properties are included in the so called Kraichnan model 
which will also be utilized in the present work. 



22j of velocity statistics. 



As opposed to a thermodynamical equilibrium, the passive scalar is maintained in a 
nonequilibrium steady state by external forcing designed to counter molecular diffusion. 
It was proved in [2^ that even in the limit of vanishing molecular diffusivity the steady 
state exists and is unique. Furthermore defining the integral scale to be infinity results in 
an infinite inertial range, divided only by the injection scale L due to the forcing. While 
the above results of the passive scalar anomalous scaling were concerned with the small 

n 

scale problem r <C L, in [15| it was observed that one obtains anomalous scaling also at 
large scales, provided the forcing is of "zero charge", qo = J (^fCiir) = 0, where Cl is the 
forcing pair correlation function. Such a forcing is concentrated around finite wavenumbers 
k ~ l/i^, which behaves similarly to a zero wavenumber concentrated forcing at small 
scales, but is more realistic for probing scales larger than the forcing scale. 



The forcing is usually taken to be statistically isotropic. Justification for this is that 
one usually expects the anisotropic effects to be lost anyway at scales much smaller that the 
for cing scale, according to a universality hypothesis by the K41 theory[r^. Nevertheless, 
in l4] it was discovered that even a small amount of anisotropy in the forcing (that can 
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never be avoided in a realistic setting) in the passive scalar equation would render the 
large scale behavior to be dominated by anisotropic zero modes responsible for another 
type of anomalous scaling. As pointed out in 14|, such behavior is nontrivial also in 
the sense that one might expect the system to obey Gibbs statistics with exponentially 



decaying correlations at large sea! 
with isotropic zero charge forcing 



es, as indeed happens for the pair correlation function 



15|. 



The purpose of the present work is to consider the small and large scale behavior of 
passive vector models stirred by an anisotropic forcing, and especially to determine if the 
phenomena of anomalous scaling and persistence of anisotopy is a general feature of passive 
advection models or just a curiosity of the passive scalar. The passive vector models arise 
as quite natural generalizations of the scalar problem and turn out to possess much richer 
phenomena already at the level of the pair correlation function. For example the pair 



correlation function of the magnetohydrodynamic equations exhibit anomalous scaling 26 1 



whereas one needs to study the fourth and higher order structure functions of the passive 



scalar to see such behavior (see e.g. 23| and references therein). It has also been argued 
that the linear passive vector models might yield the exact scaling exponents of the full 



Navier-Stokes turbulence 



a. 



The equation under study is defined as 



iii - uAui + V ■ Vui - au ■ Vvi + V^P = /j, (1) 

with a parameter a = —1, or 1, corresponding respectively to the linearized Navier-Stokes 
equations (abbreviated henceforth as LNS), the so called linear pressure model (LPM) and 
the magnetohydrodynamic (MHD) equations, z/ is a constant viscosity /difFusivity term, fi 
denotes an external stirring force, Vi is a gaussian, isotropic external velocity field defined 
by the Kraichnan model and P is the pressure, giving rise to nonlocal interactions. The 
equation was introduced in \l\ , where the authors derived and studied a zero mode equation 
for the pair correlation function in the isotropic sector and found the small scale exponents 
numerically and to a few first orders in perturbation theory (see also [4] for a more detailed 
exposition). They also reported perturbative results for higher order correlation functions 
and anisotropic sectors using the renormalization group. Although the purpose of the 
present work is to consider arbitrary values of a, some cases have already been studied 
elsewhere. The a = 1 case, corresponding to magnetohydrodynamic turbulence, has 
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probably received the most attention 



a a 



21 



25 



261]. The hnear pressure model (or 



just the passive vector modelWith a = 0. has been studied iu e.g.flflu. The liueanzed 

Navier-Stokes equation (see [17|), a = —1, was studied in [4] and numerically in [27| in 



two dimensions and is the least known of the above cases, although perhaps the most 
interesting. The above mentioned studies have been restricted to the small scale problem 
and rely heavily on the zero mode analysis, i.e. finding the homogeneous solutions to the 
pair correlation equation. For our purposes this is not enough. To capture the anomalous 
properties as discussed above, one needs to consider the amplitudes of the zero modes as 
well, as it may turn out that some amplitudes vanish. Indeed, it is exactly this sort of 
mechanism that is responsible for the anisotropy dominance in |l4| . 



We provide an exact solution of the equation for the pair correlation function with 
anisotropic forcing and study both small and large scale behavior. It turns out that for 
the "zero charge" forcing as above, the large scale behavior is anomalous even in the 
isotropic sector for all a. The anisotropy dominance seems however rather an exception 
than a rule in three dimensions, as only the trace of the correlation function for the a = 
model exhibits similar phenomena at large scales. Nevertheless, in two dimensions the 
anisotropy dominance is a more common phenomenon. Perhaps the most interesting case 
is the linearized Navier-Stokes equation for which a = —1. The field u is now considered to 
be a small perturbation to the steady turbulent state described by v. This case is unfor- 
tunately complicated by the fact that practically nothing is known of the existence of the 
steady state, although an attempt to rectify the situation is underway by the present author. 

In section [TTl we introduce the necessary tools, discuss the role of the forcing and 
present the equation for the pair correlation function in a Mellin transformed form. Details 
of it's derivation can be found in appendix |Al In section IIIII we present the solution in 
)oth isotropic and anisotropic sectors and explain the results for the passive scalar of 
l3 | in our formalism. The next three sections are concerned with the specific cases of 
magnetohydrodynamic turbulence, linear pressure model and the linearized Navier-Stokes 
equations. Although the space dimension is arbitrary (although larger than or equal to 
two), we concentrate mostly on two and three dimensions. The reasons for this are the 
considerable differences between d = 2 and d = 3 cases and the similarities of dimensions 
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d > 3. Mainly one may expect some sort of logarithmic behavior in two dimensions while in 
higher dimensions the behavior is power law like. Also the presence of anomalous scaling is 
seen to be independent of dimension for (i > 3, although the actual existence of the steady 
state may very well depend on the dimension as observed in [10]. This will be further 
studied in an undergoing investigation of the steady state existence problem. The last 
section before the conclusion attempts to shed light on the role of the parameter a as it is 
varied between —1 and 1. The actual results are collected and discussed in the conclusion. 
We also give some computational details in the appendices. 

II. PRELIMINARIES AND THE EQUATION 

All vector quantities in the equation ([1]) being divergence free results in an expression for 
the pressure after taking the divergence, 

P= (l-a){-A)-'diVjdjUi. (2) 

We may then write the equation compactly as 

Ui - uAui + Vijk {ujVk) = fi, (3) 

with a differential operator 

T>ijk = Sijdk - aSikdj + {a - l)didjdkA~^, (4) 

where is the inverse laplacian. The equal time pair correlation is defined as 

Gij{t, r) = {ui{t, X + r)uj{t, x)), (5) 

where the angular brackets denote an ensemble average with respect to the forcing and the 
velocity field. The equation for the pair correlation function is then 

dtGij — 2iy'AGij — Di^JVip^ {D^^G^p) = Gij, (6) 

where the velocity and forcing pair correlation tensors Dij and Gij will be defined below. 
The above equation should however be understood in a rather symbolic sense, as the defining 
equation for the field u is in fact a stochastic partial differential equation. The equation 
is more carefully derived in appendix |X] in Fourier variables using the rules of stochastic 



calculus. In the present form it is also very difficult to study because of the nonlocal terms 
for a 7^ 1 and the tensorial structure. We will therefore now briefly explain the structure of 
the calculations in a rather superficial but hopefully transparent way (see appendix |X] for 
details). Assuming that we have reached a steady state, i.e. dfG = 0, we rewrite eq. 
symbolically as 

- 2u'AG + MG = G, (7) 

with the effective diffusivity u' = u— ^Dm~^ and A4 is some complicated integro-differential 
operator. Taking the Fourier transform of the above equation would still leave us with an 
integral equation due to the inherent nonlocality from the pressure term. We deal with this 
now by taking also the Mellin transform (after dividing by p^), which yields 

2u'g{z) + jdz'Mz,z'g{z - z) = c{z - 2) (8) 

with a rather complicated expression for Aiz,z', see eq. flA12l) . The advantage of the above 
form is that various powers of rriy arise as poles in A4z,z' oc m^'^, the leading ones residing 
at = and z' = ^. Other poles produce positive powers of and can therefore be safely 
neglected. The residue at 2;' = cancels with the term in the effective diffusivity, leaving us 
with only the bare diffusivity u. The remaining equation can then be written in the limit of 
vanishing and u as 

- 7^ {M,,,,\z' = 9{z - = c{z - 2), (9) 

where TZ denotes the residue (the minus sign arises from the clockwise contour). The equa- 
tion is then simply solved by dividing by the residue term and using the Mellin transform 
inversion formula 

G(r-) = Jdz\r\'A,giz), (10) 

where Az is a simple z dependent function arising from the fact that we performed the 
Mellin transform on the Fourier transform of the equation. 

A. Kraichnan model 

We define the Kraichnan model as in with the velocity correlation function 

=: 5{t)D,,{r-m,) (11) 



where we have defined the incompressibihty tensor Pij{q) = Sij — qiQj and denoted (f^g := 
(0- Defining 

DmM)= , ,J%2^,/2 (12) 
and applying the MeUin transform (See e.g. 13| and the appendix A of [23]) we have 

oo 

51 (q) := / '^^^'^'^rn. = dm. {z')q-^'~\ (13) 



where 

and z' is constrained inside the strip of analyticity ~d < Re{z') < ^. The parameter ^ takes 
values between zero and two and measures the spatial "roughness" of the velocity statistics. 
We observe that the scaling behavior of the correlation function is completely encoded in 
the pole structure of Mellin transform, with e.g. the pole at 2;' = ^ corresponding to the 
leading scaling behavior of the velocity structure function. 



B. Decomposition in basis tensor functions 

Being a rank two tensor field, the pair correlation function may be decomposed in hy- 
perspherical basis tensor functions as in . Such a decomposition is also an important 
tool in analyzing the data from numerical simulations, as witnessed e.g. in jj]. We shall 
be concerned only with the axial anisotropy, and apply this decomposition on the Fourier 
transform of the pair correlation function. This has the advantage of making the incom- 
pressibihty condition very easy to solve, among other things. We consider only the case of 
even parity and symmetry in indices, which leaves us with a basis of four tensors: 



' 4(p) 


= \p\ 






= \p\' 


~^didj^\p) 




= i^r 


\pidj+pjdi)^\p 




= i^r 





(15) 



with the actual decomposition 

G,,{p):=Y,B^i{p)G\{p). (16) 

b,i 
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Here is defined as := where is the hyper spherical harmonic function 

(with the multi-index m = 0). It satisfies the properties 

A^\p) =0 

p.V^\p) =l^\p). (17) 
The same decomposition will naturally be applied to the forcing correlation function as well. 

C. The forcing correlation function 

We require the forcing correlation function to decay faster than a power law for large 
momenta and to behave as Cij{p) oc L'^{Lp)'^^ for small momenta with positive integer 
A^. The = case corresponds to the usual large scale forcing with a nonzero "charge" 



go = JdrCi 
large scales 



r) and is responsible for the canonical scaling behavior of the passive scalar at 



151 ] , whereas any > corresponds to a vanishing charge 151] • Applying the 



Mellin transform to such a tensor (decomposed as above) yields 

oo 

C!,{p) =y"^^^+^a,(n;p) = ]p]-'^-^5^4(p)cf(z), (18) 





with 

c^(z) = ' , Ueiz) >-d~ 2N, (19) 

and the strip of analyticity —d — 2N < TZe{z). The details of the actual cutoff function 
are absorbed in the constants and play no role in the leading scaling behavior. All the 
interesting phenomena can be classified by using only the cases A^ = and A^ = 1. We will 
mostly be concerned with the latter type of forcing which is also of the type considered in 
Q, 3- By inverting the Mellin transform we would obtain an expression for the forcing 
correlation function 

a,{t,r) =ldz\r\^c^{z)K^\z)Bl^{r) (20) 

where the matrix K is defined in appendix [Dl We note that determines the large scaling 
behavior of the above quantity as r"'^ or r"^"*^, depending on the forcing, while the matrix K 
is responsible for the small scale behavior oc r', where / is the angular momentum variable. 
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D. Mellin transformed equation and overview of calculations 



As mentioned earlier in this section, equation ([6]) is much too unwieldy for actual com- 
putations. In appendix |X] we perform a more careful derivation of the equation in Fourier 
variables and by using the Ito formula. The resulting equation (]A9p still has an inconvenient 
convolution integral. By applying the Mellin transform, we obtain an equation 

-Dm-%{z) - Do\gb{z - + j^z' draSz')^%z' 4+z~z'9c{z - z') 

= Cb{z - 2). (21) 

for the Mellin transformed coefficients cjb of the tensor decomposition f|T6l) (defined explicitly 
in eq. (lAlip ). The matrix T is defined in eq. (1A13P and involves rather difficult but 



manageable integrals, and A is defined in eq. f lA7p . The integration contour with respect to 
z' lies inside the strip of analyticity TZe{z) < TZe{z') < 0, determined from eq. (lA14p . For 
small values of m„ the contour may (and must) be completed from the right. The reason for 
performing the Mellin transform becomes evident when one studies the pole structure of the 
functions dm^{z') and T: first two (positive) poles occur at z' = (from T) and at 2;' = ^ 
(from dm^{z')) and correspond to a term oc and a constant in m^, respectively. The 
former of these cancels out from the equation, hence one is free to take the limit 0. 
This leaves us with a simple equation 

-XUz -0- ^dH4+z^f9ciz - = ^Cbiz - 2). (22) 

From now on we absorb Dq in the functions Cb- In appendix [B] we have applied the in- 
compressibility condition to the correlation function Gij{p) and the equation, which has the 
effect of leaving us only with two independent functions to be solved, gi and g2. Applying 
also a translation z ^ ^ + ^ in eq. fl23l) , we have 



with the definitions 



(Al + A + B-X)h(2) = f(^ + e-2), (23) 



and 



-d+^,d+z 



h 


= (^1,^2)"^ 


f 


= (Cl,C2)^ 
















(24) 



(25) 



All that remains now is to invert the matrix equation, although in the isotropic sector and 
in two dimensions it reduces to a scalar equation. 



III. THE SOLUTION 



Inverting the Mellin and Fourier transforms enables us to write the full solution as 

G.,(r) = -jdz\r\'h^{z)P^K ■ B,,{r), (26) 

where we now have a projected version of the matrix K due to the incompressibility condition 
(see appendix [D|) . and 

h(z) = - (^Al + A + B-x) ^f(2 + e-2). (27) 

The strip of analyticity is now 

2-d-^-2N< ne{z) < 0, (28) 

where iV = for the traditional nonzero charge forcing and = 1 for the zero charge 
forcing. We should note that there may in fact be poles inside the strip of analyticity due 
to the solution h, which is just a reflection of one's choice of boundary conditions. 



A. Isotropic sector 

In the isotropic case when I = 0, we have Blj = 6ij, Bfj = riVj and the other B's are 
zero. The equation of motion (l23l) is now a scalar equation, hence we only need the (1, 1) 
-component of the matrix. 



(^Al + A + B ■ x) = -{a - l){a^ - 1 - a - d)T {1 + i/2) T (1 + d/2) 

where the equality applies up to a constant term that will be absorbed in the forcing, and 
we have defined the polynomial 

Pa{z) = {a-lf{l + d)iz 

+ {z + i- 2) {{d - l)z{d + z) + a{a{d - l)d + 2z) ^ + a^{d - (30) 
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This is the same expression (only in a sUghtly different form) as in The expression (126|) 
for the inhomogeneous part of the correlation function becomes 

G.,(r) =jdz\rY^a{z)c,{z + ^ - mM)^^^. (31) 

where we have introduced the incompressibility tensor 

V^J{Z) = [{z + d- l)5ij - ZViVj] (32) 

and irrelevant constant terms were absorbed in the forcing Ci. Henceforth such an assumption 
will always be implied unless stated otherwise. 



B. Anisotropic sectors 

Now the task is to find the poles of the inverse matrix of ^Al + A + B ■ , that are 
completely determined by the zeros of its determinant. Denoting 

M := A + B . X = - - - + - '^'^A , (33) 

^ + ^ \ri2 - r42 - (/ - \)r^2 + (l - l)ln2 J 

where r and A are defined in appendix O, we may write 

det (Xl + m) = A^ + A trM + det M. (34) 

We refrain from explicitly writing down the determinant, since the full expression is rather 
cumbersome and not very illuminating. It may however be easily reproduced by using the 
components Tij given in appendix O 

C. Two dimensions 

The two dimensional case deserves some special attention. From the incompressibility 



requirement in eq. fIBll) and by direct computation using the two dimensional spherical 
harmonics oc e'^, one can see that the correlation function satisfies the propotionality 

G,,ip)^{f-l{l-l)f)P,,ip). (35) 

Therefore in two dimensions the equation is a scalar one also in the anisotropic sectors. A 
formula for the solution then becomes 

G„(r) = -|&|rrp^^^|^(P-K)-4(f). (36) 
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where F = Al + A + B ■ X. 



D. Example: Passive Scalar 

As one of the main themes of the present work is to consider the effect of a forcing locahzed 
around some finite wavenumber ruf (x 1/L instead of zero, it is useful to review the case in 
[l^ by the present method (see e.g. 12, 19, 23] for more on the passive scalar problem), even 
more so as the magnetohydrodynamic case in two dimensions bears close resemblance to the 
passive scalar (indeed the two dimensional case can be completely described as a passive 
scalar problem with the stream function taking place of the scalar). Using the methods 
above, we arrive at an expression similar to (12^ . 

G(r) = Y: J"^ l-l' °"''Ml""/( -^S l y (37) 

where we have again written the generic constant C in which we will absorb finite constants. 
In the above equation, N equals zero or one corresponding to the nonzero and zero charge 
forcings and 

Mz) = {d - - z){l + z + d + ^ - 2) + e/(/ - 1). (38) 

The strip of analyticity is now —d — ^ < TZe{z) < 0. Consider now the isotropic sector I = 
with the nonzero charge forcing, i.e. N = 0. We then have (neglecting the zero modes) 

G«(x)^C'L-^/f.|r/i|-^jjML. (39) 

For r <^ L the integration contour must be completed from the right, thus capturing the 
poles z = 0, z = 2 — C,, ■ ■ ■■ The small scale leading order behavior is therefore 

+ (40) 

where the dots refer to higher order powers of r. The large scales r ^ L require a left hand 
contour, resulting in another scaling regime, 

Gz=o = C'^^^^L'r'-'-' + ... (41) 

We note that the above solution is constant at r = and zero at r = oo, thus satisfying 
the boundary conditions. We conclude that the solution is completely nonanomalous, i.e. 
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respecting the canonical scaling. 



Consider now instead the zero charge forcing with = 1 that is localized around 
p = 1/ L instead of p = 0. The large scale pole due to the forcing at 2; = —d — ,^ + 2 cancels 
out and we are left with 

G,=o = C'jdz \r/L\' -r ■ (42) 

There is now no large scale scaling behavior (the decay is faster than a power law). By 
looking at the / = 2 sector, 

(with a different generic constant C"), we see that the relevant scaling behaviors are obtained 
from a solution of the equation 

^Ij2{z) = d\-2 + z)- z{-2 + z + i) + d{-2 + z){~l + z + = 0, (44) 

giving the large scale behaviour of the / = 2 sector with the exponent 



(45) 

Therefore we conclude that the large scale behavior is dominated by the anisotropic modes. 
Note that the anisotropic modes are also anomalous in that they are not obtainable by 
dimensional analysis. 

IV. MAGNETOHYDRODYNAMIC TURBULENCE 



Settirig a = 1 in eq. ([I]) yields the equations of magnetohydrodynamic turbulence (see 



e.g 
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26j and references therein) . This is a special case in that the problem is completely 
local due to the vanishing of the pressure term. In practical terms, the quantity A is zero, 
hence we only need to consider the zeros of the determinant of M in eq. flMl) . 

A. Isotropic sector 

The isotropic part of the correlation function becomes 

G.,(r) = Cy^z|rr^ii^i^P.,^ (46) 
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where 



Po{z) = {d- l)z{d + z) + {{d - l)d + 2z)^ + {d- 1)^2 (47) 
with another generic constant C. We find the usual poles at 
Zn = 2 - ^ + 2n, 



where n is a nonnegative integer. For the nonzero charge type forcing we have c1{z + C,—2) oc 
l/{z + d + ^ — 2), which presents another pole. On the other hand, for the zero charge forcing 
we have c\^{z + — 2) oc l/{z + d + which cancels with a zero of the gamma function. 
It turns out that this sort of cancelation occurs for each model, rendering the large scale 
behavior anomalous. We will postpone the arbitrary dimensional case until the end of the 
present sector and instead consider first the three and two dimensional cases. 

B. Anisotropic sectors 

Note that since detM oc ^'i^d+z,d+^^ ^^e inverse of M is only proportional to Xj'^^^^^^^, 
so the correct form to look at is actually det M/ Xi^d+z,d+^- Dropping z -independent terms 
we have 

detM r(^)r( '+-+f^-^ ) ^ ^ 



Xi+d+.4+^ r (tt^) r (^^^) 

where "^liz) is a fourth order polynomial in z and C is a z -independent constant. Due to its 
rather lengthy expression, we shall consider the whole problem in three and two dimensions 
only. We note immediately that there's also an infinite number of solutions due to one of 
the gamma functions, namely at 

z = l + 2-^ + 2n (50) 

for nonnegative integers n and even /. The other gamma function cancels with the terms 
from flD2l) . 



14 



(a) d=3 (b) d=2 

FIG. 1: The MHD scaling exponents of the isotropic, I = 2, I = 4 and 1 = 6. In (a) The isotropic 
poles z+ > z— are adjoined at ^ = 1. The dashed line in (a) corresponds to the forcing with 
nonzero charge with a pole at —1 — ^, whereas for the zero charge forcing there are no poles. In 
(b) the zero modes are never adjoined. 

C. d=3 

We have the four solutions to "^1(2) = of which the following two are dominant in the 
small and large scales, 



^r = -^^±yAT2i2-OVB, (51) 



where 



A =(2 + 0(2/(/ + l)-6-0 + 17, 

B =(2 + 0(2/(/ + l)+0 + l, (52) 



which match exactly to the results obtained in [6|, |25[ , after some convenient simplifications. 
The isotropic zero modes are 

^± = ^(3-^± V3(l-0(3 + 0)- (53) 

We have plotted the leading poles in Fig. ([T]) from / = to / = 6 together with the pole due 
to the nonzero charge. We note that the isotropic exponents become complex valued for 
^ > 1, implying an oscillating behavior and therefore a positive Lyapunov exponent for the 
time evolution [lO|, l26| . The above steady state assumption therefore applies for < ^ < 1 
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only in the isotropic sector. The fact that the anisotropic exponents are continuous curves 
for all < < 2 seems to imply that the steady state exists for all ^ in the anisotropic 
sectors. Indeed, in ^ this was shown to be the case by preforming a more careful eigenvalue 
analysis. 

1. Nonzero forcing charge 

In the isotropic sector for the forcing with nonzero charge = we have 

G,j[r)\i=o =-CL ^ dz- V ^1 ^nrP+^+g^ ^^^> 

J [z - z^)[z - z+)[z + I + i)V [^^) 

with the contour bound — 1 — ^ < lZe{z) < and z_ < —\ — ^ < Zj^ < C again denotes 
some generic finite (and positive) constant. The pole 2;+ divides the strip of analyticity in two 
parts, which correspond to different boundary conditions. Small scale behavior corresponds 
to picking up the poles to the right of the contour and large scale behavior corresponds to 
left hand side poles. We note that both the zero modes z±. are negative, except that 2+ = 
at ^ = 0. Therefore 2;+ cannot be a large scale exponent, as the solution has to decay at 
infinity. The real strip of analyticity is then in fact — 1 — ^ < TZe{z) < z^, thus resulting in 
the small scale behavior 

G< = C,r^^n,{z^) (55) 



and the large scale behavior 



G>. = C2r-i-%(-l-0 (56) 



We note that the large scale behavior is determined by the forcing and therefore respects 
canonical scaling. 

2. Zero charge forcing 

Because of the pole cancelation we now have a similar expression, 

G.,(r)b.o =C'L^-^hz '"^'^'J^^'^'^'^^Pji. (57) 

J [z - z_){z - z+)r {^j^} 



16 



with the strip of analyticity is now —3 — ^ < TZe{z) < 0. The contour bound now encloses 
both the zero modes (see again Fig. ([1])). In addition to the above considerations with a 
forcing of nonzero charge, we conclude that cannot be present at small scales due to 
regularity conditions at ^ = 0, so the real strip of analyticity is in fact z_ < TZe{z) < z+. 
This gives rise to the small scale behavior 



G< = Cir^+P.,(^+) (58) 



for the small scales and 



G> = C,r^-V,,{z.) (59) 

for the large scales. The large scales are therefore dominated by the smaller zero mode 
Z- instead of the exponent — 1 — ^ as with the nonzero charge forcing and is therefore 
anomalous. However, unlike in the passive scalar case, the anisotropic exponents are sub- 
dominant at both small and large scales (see Fig. ([1])) and we therefore conclude that there 
is isotropization at both scales. 

D. d=2 

The (dominant) zero modes in two dimensions are 



zt =-4-e+ V4/2(l + 0+e 

=_2,i-^Al\l + 0+e (60) 



of which we separately mention the isotropic zero modes, 

z^ =-2- e (61) 
The expression for the inhomogeneous part of the correlation function is 

G.,(r)|,=o =C'L'''^jdz\T/LYcL{z + i-2)V.,{z)^^^^ (62) 
with contour bound is —2 — ^ < 2; < together with the bound from the forcing. 
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1. Nonzero charge forcing 



Because cl oc 1/(z + ^), the expression for the isotropic sector of the correlation function 
simphfies to 

G.,(r)|,=o =C'L'-^jdz\r/LYV.,{z)^-^^^^ (63) 

where the bound is now — ^ < 7le{z) < 0. Note the appearance of a double pole at 2; = — ^ 
giving rise to logarithmic behavior. There are now no poles inside the contour bound, so 
finding the asymptotics is easy. We observe that there are no small scale poles and therefore 
the correlation function decays faster than any power at small scales, whereas at large scales 
we have 

G>. = C"log(r/L)L2-«|r/Lr%(-0 + C L'-^\r / L\-^V[^{-0, (64) 

where V[j{—^) = 6ij — TiVj to ensure incompressibility and other next to leading order 
nonlogarithmic terms were discarded. By looking at Fig. ([1]) we see that there is a hierarchy 
of small scale exponents in the anisotropic sectors. We therefore make the conclusion that in 
two dimensions the anisotropic effects in the MHD model are dominant at small scales for a 
forcing of nonvanishing charge, conversely to the passive scalar case. Note that setting = 
in the above equation reproduces correctly the usual logarithmic behavior of the diffusion 
equation steady state with an infrared finite large scale forcing. 



2. Zero charge forcing 

We now have oc l/(z + ^ + 2) and the isotropic correlation function becomes 

G.,(r)|,=o =-C'L'-^ldz\r/L\^V.,iz)^^^^ (65) 

with the usual strip —2 — < TZe{z) < 0. There are no double poles and the leading simple 
poles are just at z = —C, and z = —2 — C,, so the asymptotic behaviours at small and large 
scales are simply 

G<- = C'|r/L|-%(-0 

Gf, =C'\r/Lr-^V.,{-2-0. (66) 
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As in the three dimensional case, all the anisotropic exponents are now subleading at both 
small and large scales (see Fig. ([1])), so we conclude that there is again isotropization at 
both regimes. Note also that the large scale behavior is due to the forcing and therefore 
nonanomalous. 



3. Any dimension, zero charge forcing 

For the sake of completeness, we write explicitly the solutions in any dimension d > 2 in 
the isotropic sector for the zero charge forcing: 



2 z+ - I 2+rf+^++g 



(2-^-^_)(2-e-^+)r(2 + d/2)^ ^ ^ 



2 z+ - z_ Y ( 2±^±I^ 



where we have neglected the possible exponentially decaying terms. The anisotropic sectors 
produce rather cumbersome expressions and we will be satisfied with only the numerical 
results in the figures. We observe that the large scale behavior is always dominated by the 
negative zero mode exponent z^ and is therefore always anomalous (except in two dimen- 
sions). It is also fairly easy to see that the anisotropic exponents are always subdominant, 
so that there is isotropization at both small and large scales. 



V. LINEAR PRESSURE MODEL 

Setting a = in eq. ([1]) produces the equation known as the Linear Pressure Model 
(LPM) (see e.g. and references therein; sometimes this model is just called the 

passive vector model) By looking at equation (1A9I) . we see that when G oc d'^^^p), the left 
hand side evaluates to oc a^|pP~'^~'>P,,, (p). Therefore for a = there is a constant zero mode 

nn 

analogously to the passive scalar case. This is true for the anisotropic sectors as well ([2, IJ])- 
This constant zero mode however vanishes for the structure function, so in the present case 
we also consider the next to leading order term. The first thing to note in the isotropic 
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sector is that when a = 0, z = —d is a solution of the equation 

1 



;(rf + i)r(e/2)r(i + c//2) 



r(-./2)r(^)r(i±f^) _ 
+^o(.) — r(^)r(-f^) °' ^^^^ 

where 

Po(^) = ^(^^'(^ + e - 2) - ^(^ + e - 2) + - 2)(^ + e - 1) - 0- (69) 

However, as we see from the definition of the incompressibihty tensor in eq. (152]) . for the 
trace (in indices) we have 

V,,{z) = {d-l){d + z), (70) 

which produces a cancehng z + d term in the numerator. A physically more realistic quantity 
would however be a contraction with x'^x^ than the trace, since we are more interested in the 
structure functions of the model. Another exact solution is 2; = 2 — ^. Other nonperturbative 
solutions can only be obtained numerically. 



A. Any dimension 

We have plotted some of the poles in Fig. ([2]) in three dimensions. Remembering the 
z = —3 solution, we see that the anisotropic exponents are less dominant with increasing / 
(a behavior repeated for higher / as well). 



1. Nonzero charge forcing 

The contour bound is now 2 — d — ^ < TZe{z) < 0, so there is no controversy in the choice 
of which poles to include. The small and large scale behaviors are similar to the passive 
scalar, and for completeness, we give the results in any dimension: 

G< = AL2-%-Sr2-%(2-0 

G'>. =A'L'-^\r/L\^~''-^Vij{2-d-0- (71) 
The A, B and A' are somewhat complicated transcendental functions of d and ^. 
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(a) d=3 (b) d=2 

FIG. 2: The Linear Pressure Model scaling exponents of the sectors Z = 0, / = 2 and / = 4 in 
three (a) and two (b) dimensions. The z = and z = —d are omitted for the sake of clarity. We 
note that in two dimensions, there is a z = 2 exponent in the / = 2 sector but the / = 4 sector's 
exponent goes slightly above z = 2. 

2. Zero charge forcing 

Now the forcing contributes a pole oc l/(z + ^ + (i) and the contour bound is —d — ^ < 
TZe{z) < 0. The quantity 70 in eq. (1^5]) has a zero there that cancels with the pole from 
the forcing. Therefore we again conclude that the forcing doesn't contribute in the scaling. 
The small scale behavior is therefore same as above, but the large scale isotropic sector of 
the correlation function behaves as 

G> = C'\L\'-^ {A'\r/L\-''V,,{-d) + B'\r/LrV,,{z_)), (72) 

where A' and B' are again some nonzero constants (depending of d and ^), Z- is the / = 2 
large scale mode (see Fig. ([2])) and we have the traceless tensor 

Vij{-d) = dviVj - Sij. (73) 

By looking at Fig. we observe that the subleading exponent z_ is smaller than the 
anisotropic exponent Z = 2 in three dimensions and Z = 4 at two dimensions (except when ^ 
is close to two, when the / = 2 exponent is larger than the / = 4 exponent). Therefore the 
trace of the correlation function is dominated by the anisotropic modes. 
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(a) d=3 (b) d=2 

FIG. 3: The Linearized Navier-Stokes equation exponents for sectors I = 0, I = 2 and / = 4 (the 
legend apphes to both figures) at three and two dimensions. In (a) the I = 4 curves run shghtly 
below and above the curves z = —3 — and z = 2, respectively. Other than leading exponents are 
also displayed. 



VI. LINEARIZED NAVIER-STOKES EQUATION 



Setting a = — 1 in eq. ([T]) yields the Linearized Navier-Stokes equation (see e 



•g. HQ, 



24l|). The equation may be considered as zeroth order perturbation theory of the full Navier- 
Stokes turbulence problem, from which one can at least in principle proceed to higher orders 
in perturbation theory. It will also serve as a stability problem where the background flow is 
determined by the Kraichnan ensemble instead of a solution to the Navier-Stokes equation 



m 



see chapter III of 2J]). Not much is known of this case, except for the perturbative results 



^. The eq. fl29|) becomes 

1 



d(d + 2 + or(i + e/2)r(i + (i/2) 



'J ^ r(-./2)r(^)r(^) 

with 

p_i(z) ={z + ^-2) {d\z + - + 0' + d{z - l)z + d(e - 1)0 - 4(rf + 1)^^. (74) 

We choose to save space by not writing down explicitly the determinant for the anisotropic 
sectors. The expression may be reproduced by using the results of appendix O We will also 
refrain from explicitly writing down expressions for the correlation functions, as it turns out 
that whichever sector has the leading exponents varies quite a bit with different values of ^. 
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A. d=3 with zero charge forcing 



The contour bound is, as usual, —3 — ^ < TZe{z) < and again one observes a cancelation 
of the corresponding pole. Inspecting Fig. one observes quite wild behavior of the various 
scaling exponents at a first few sectors. A notable similarity to the three dimensional MHD 
case (a = 1) are the exponents starting at and —3 and joining at ^ ~ 0, 35. However in 
the LNS case one also sees similar behavior near ^ = 2. Indeed one is tempted to assume 
the existence of a steady state only for ^ near zero and two. The same conclusion could 
be drawn for the anisotropic sectors as well. We will further discuss this at the end of the 
paper. We will be satisfied with only reporting the scaling behaviors as the procedure for 
finding them is close to above cases. Assuming the steady state exists for ^ close enough 
to zero and two, we conclude that for ^ near zero, the small and large scale are dominated 
by the isotropic exponents starting at and —3, respectively. For ^ near 2, one instead 
observes I = 2 dominance at small scales and I = 4 dominance at large scales. We have 
deliberately neglected the nonzero charge forcing, as that would only bring about the familiar 
nonanomalous — 1 — ^ scaling at large scales. 



B. d = 2 with zero charge forcing 

The behavior of the scaling exponents are much nicer, as can be seen by looking at Fig. 
([3]). For < ^ < 1,3, we see the small scales dominated by the / = 4 anisotropic sector, 
and the large scale by the / = 2 sector. For other values of ^ the / = 2 anisotropic sector 
dominates the small scales as well. The I > 4 anisotropic exponents are all subleadingwith 
respect to the ones in the figure, and indeed respect the usual hierarchy of exponents p|. In 
any case, the isotropic exponent is subleading. 



VII. THE EFFECT OF VARYING THE PARAMETER a 

It is useful to discuss also other values of a beside the discrete values a = 1, 0, —1. More 
specifically, looking at Fig. (jl]) we see how the closed contour determining the leading scaling 
exponents is deformed as a varies from a = 1 and a = —1 to 0. Both end up as curves z = 
and z = —3. Also, as we know that when a = 1 the steady state exists for f < 1 in the 

n 

isotropic sector [26i] (and for all ^ in the anisotropic sectors [6,]), it now seems even more 

23 



(a) (b) 

FIG. 4: The leading isotropic exponents as a is varied from 1 to (a) and from —1 to (b) in 
three dimensions. The darkest curves correspond to a = 1 and a = —1. 



reasonable to expect the steady state to exist for all ^ in the a = case. 

VIII. CONCLUSION 

The purpose of the present paper was to present an exact solution for the two point 
function of the so-called a -model in the small and large scaling regimes, which incorporates 
the magnetohydrodynamic equations, the linear pressure model and the linearized Navier- 
Stokes equations. The phenomena of anomalous scaling and anisotropy dominance was 
investigated in each model with emphasis placed in the zero charge forcing concentrated at 



a finite wavenumber ~ 1/L as in 15| . Below we briefly summarize the findings in each model. 



For the magnetohydrodynamic equations with a = 1 the leading scaling behavior 
was observed to be anomalous and isotropic at both small and large scales in three 
dimensions for the zero charge forcing, in accordance with previous small scale results 

0,fl~ 



25| . In two dimensions with nonzero charge forcing one observes anomalous and 
anisotropic behavior at small scales, while the large scales are dominated by logarithmic 
behavior. The mechanism of the small scale anisotropy dominance is strikingly similar 
to the passive scalar large scale anisotropy dominance, except that in the MHD case the 
phenomenon results from the nonzero charge forcing. The zero charge forcing case in two 



dimensions is in agreement with the results in 26 1. 
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For the linear pressure model with a = and zero charge we recovered the small 
scale exponents of [2]. The small scale behavior is now dominated by the isotropic and 
canonical scaling exponent z = 2 — ^ (neglecting the constant mode by considering the 
structure function). The large scale behavior was seen to be dominated by a curious 
isotropic zero mode z = —d, although the trace of the structure function exhibits anomalous 
and anisotropic behavior at large scales. The nonzero charge forcing simply renders the 
large scale behavior canonical. The existence of the steady state is nevertheless controversial 
in two dimensions and requires further study. 



The linearized Navier-Stokes equations corresponding to a = —1 seem to be the 
most interesting of the models considered, even more so because it is also the least well 
known. There still remains the question of the existence of the steady state, without 
which one cannot claim to have completely solved the problem. One may however 
conjecture it's existence at least for small enough C, (at least in the isotropic sector), 
in which case the small and large scales are dominated by the isotropic anomalous 
scaling exponents in three dimensions. In two dimensions, the small scale exponents 
coincide with the somewhat rough numerical estimates of 27|, the difference now be- 
ing the absence of the scaling oc due to the forcing. Indeed, it was observed that 
both the small and large scales were dominated by anomalous anisotropic scaling exponents. 



Although the linear equations above with the somewhat crude Kraichnan model are 
certainly some distance from the real problem of turbulence, similar scaling behavior has 
been observed in real and numerical simulations (see e.g. 0, Q] and references therein), 
namely implying that the scaling exponents in each anisotropic sector are universal as 
outlined above. Probably the closest case to the real Navier-Stokes turbulence is the 
linearized Navier-Stokes equation. The equation arises usually as one tries to verify the 
stability of a given stationary flow by decomposing the velocity field as v + u, where v is 



the stationary, time independent term and u is a small perturbation [2J]. If one can show 
that u decays in time, the velocity field v is indeed a laminar, stable flow. In our case 
V is determined by the Kraichnan model and we are now concerned with the stability of 
the statistical steady state. It has been pointed out in [4] that in such a case one might 
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be able to show that higher order perturbative terms are irrelevant in the sense of the 
renormalization group, thus implying that the steady state is in fact in the same universality 
class as the full NS turbulence. This would mean that the anomalous scaling exponent of 
the linear model is equal to the NS turbulence exponent. All this would of course depend 
on the existence of the steady state for u. Unfortunately it seems that such a steady state 
does not exist for the exponent ^ = 2/3, which could be a sign of incompleteness of the 
Kraichnan model or a symptom of the general complexity of the problem of turbulence. 
The stability and existence problem will be studied more carefully in a future paper by the 
present author. 
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APPENDIX A: EQUATION OF MOTION FOR THE PAIR CORRELATION 
FUNCTION 

We take the Fourier transform of equation ([1]) and rewrite it as a stochastic partial 
differential equation of Stratonovich type as 

duiip) = up'u,{p)dt - V^^,jd^qdV,iq) o u^{p - q) + dF,{p), (Al) 

where we have dropped the t -dependence and denoted 

= ^ i^iaPb - aSibPa) + i{a - l)piPaPb/P^ (A2) 

and defined the Stratonovich product 

dVy{q)ou^{p- q) = dVy{t,q)%{t + ^,p ~ q). (A3) 



As argued in 281] by physical grounds, the symmetric prescription ^(0) = 1/2, corresponding 



to the Stratonovich definition of the SPDE, is the correct way of defining the equation. We 
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will however use the relation it^(t + ^,p) — u^{t,p) + \du,j,{t,p) to transform the equation 
into a following Ito SPDE, 

du,{p) = vp''u,{p)dt+]^%,jd:'q%M)K-^Mp) 

... - Vl,j(fqdVM)u,{p -q) + dF,{p), (A4) 

where we have used the relation 

dVi{t,p)dVj{t,p') = %{,p)5\p + p')dt. (A5) 

The first integral on the right hand side of the Ito SPDE can be done exphcitly, resulting in 

j<^<lD^M)K~>M = Dm-^P'UP) + h'-%{p) + (A6) 
where the incompressibility condition piUi{p) = was used, and denoting 
A = (a - 1) (d + 1 + a(l - 0) ^<^^f</2)r(y2)c. 

Applying the Ito formula to the quantity 

{u,{t,p)u,{t,p')) = G,,{t,p)5\p + p') (A8) 

and by assuming stationarity, one obtains the nonlocal PDE (with obvious p dependence 
omitted) 

[2v - Dm-^] IpI^G,, - XDo\pt^G,, + V^^VJ^, jd'qD,M)G,p{p - q) = Q,- (A9) 
Using the SO(d) decomposition for G, 

(AlO) 

a 

(and similarly for C), dividing the equation by p^, and by taking the Mellin transform of 
the equation while remembering the definition 

oo 

Gm = I ^w'+^G,j{wp) = IpI"'^-^ Y1 BUP)9aiz), (All) 

^ 1 

and by expressing D in the integrand as an inverse Mellin transform, we finally obtain the 
equation 

[2u - Dm-^] g,iz)-XDoU^ - + jdz'dmS^')i:%^,^,^^_^,Uz - z') 

^c,{z-2), (A12) 
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where we have defined (note the transpose in definition) 

P,^ip-q)Bliq) 

b 

with the strips of analyticity, 



^d+z'^+z-z'-l^ijKP) - \P\ ^itxv^jpa /" .d+z' I \d+z-z' ' K^^'^) 

J \p-q\ \q\ 



ne{z) - ne{z') < 
ne{z') < 

d + 7^e(2) > (A14) 

such that the 9x9 matrix T is independent of p. The matrix elements T*'^ can be determined 
exactly by computing the right hand side integral, which is the subject of the next appendix. 
As mentioned in sec. IIIDl the first poles on the right occur at = and z' = ^, which 
results in the equation in the limit of vanishing m^: 

2ugt{z) - DoXUz - + d^Sm'^Uz) - D,T%^^^^,^^{z)Uz - 

= Cb{z-2). (A15) 

We have defined the residue matrix 

R'"^ = 7^ {T%,^d+z-z) l.'=o (A16) 
and used the residue of the velocity correlation at z' = ^: 

7^.'=^ = -^0. (A17) 

APPENDIX B: INCOMPRESSIBILITY CONDITION 

The incompressibility condition for u and / amounts to requiring that the contraction of 
the covariances f|T6|) with p is zero, i.e. 

= {pfgi + ipfgz + PfgA) ^\p) + IpI' {{i - 1)928^ + g.d,) ^\p) 

= 0, (Bl) 
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which gives a system of equations 

9i + I93 + 94 =0 

{l-l)-g,+g, =0. (B2) 

We can achieve this conveniently by defining a projection operator 

1 0^ 




(B3) 



where 













IX. f) 









X= ) \. (B4) 

-1 /(/-I) 

The solution to eq. flB2l) (and a similar one for the forcing) can then be written conveniently 

as 

I u \ / f \ 

(B5) 

We also rewrite the matrices R and T in block form as 

f Ri R2\ /^A b\ 

R = ; T,+5,,+, = . (B6 

\^R3R4y \CBJ 

Note the above definition of T with a translation z z + ^. R is independent of z. By 
operating with P on eq. (123!) . we obtain the equations (after translation z ^ z + 

[2u-Dm-^] h{z + 0+dmM (Ri+ R2-X)h(^ + ^ - ADoh(^) - Do (A + B ■ X) h(^) 
= f(z + e-2) (B7) 

and an identical one but multiplied by X from the left. Thus we see that we only need the 
upper 2 by 2 matrices from T. By using the definition eq. (1A16|) and the results for T°* in 
appendix O we obtain 

which results in a cancellation of the remaining mass dependent terms. The remaining 
equation depends now only on the physical diffusivity z/. Solving the equation iteratively 
would amount to a series expansion in powers of u or z/"^, but we shall only consider the 
1/ — i> limit, which produces the solution in eq. (|23|1 . 
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APPENDIX C: NECESSARY COMPONENTS OF THE MATRIX T 



Due to incompressibility, only some of the components of T will be needed. Computing 
the integrals of the type in flAlSp can be performed by using the result 

A , f =: A.«,., \pt''-^'' ■ e), (CI) 

where we have denoted hy q-e the angle between q and the z -axis and defined 

_ T{d/2 + / - a)V{d/2 - f3)V{a + (3- d/2) 

The tensorial structure can be obtained by partial integrations and by taking derivatives in 
p. We will further define (note the transpose in the definition) 



rpafe _ ^l+d+z,d+£ab/\ fn'i\ 



The necessary components of r are (others do not contribute due to the incompressibility 
condition) : 



-11 



[l + a^){d-l){l- z) -a^^{z + d + ^-l) 



{I- z-Oil + z + d + ^-2) 



T 
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{l-z-Oil + z + d + ^-2) 



.21 



.22 



1 + z + d-2 



^ 'l + z + d + i 
{d-l){l + z + d-2) 



+ 



d-l+i- 
2)aa\l 



-l+d+i+2 

z-l + 2 
3) + 2a{z + d + l) + 1 -3) 



1 + z + d + ^-2 ' {I- z 

(g- 1)^(2-0^(^^-5/ + 6) 
{l-z-2){l + z + d + i-2){l + z + d + i-A) 



2){l + z + d + i-2) 



2al^{z + d + i-l) 2 f z + l{d + i-l)-{d + i){z + i) 

{l-z-i){l + z + d + i-2)^ " ^ 

(/-i)e(^+e-i) (/-i)(2-o(rf+oe \ 

{l-z-C){l + z + d + i-2) {I - z - 2){l - z - + z + d + i - 2) ) 

a^d-l + a{l-2) + z + i)-d-l 4a(/-2)(2-Oe 

^ {I- z-i){l + z + d + i-2) {l-z-2){l-z-i){l + z + d + i-2) 

2e(a-l)2(/2_5/ + 6) (2-0 
{I- z-2){l- z-i){l + z + d + i-A){l + z + d + i-2) 
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41 amd + ^){z + 0-l{d + ^-l)-z)-2a^ 

= J c ^ 

I — z — t, 

^ ^ (i + 1 + 2a{l - 1) + a\l + 2d-cP- i{d - 1)) ^ {a-lf{d + l)(2-i)i 



-(2 -Of 



^42 



a - 1)2 (d2 + /(/ + !)) + (/ ((1 + 2/)(l - 2a) + a2(l + 3/ - f)) 
{I- z-2){l + z + d){l- z-^){l + z + d + ^-2) 

2(a - l)\d + l){d+l + l)i{e -6^ + 8) 

' {l-z-2){l + z + d){l-z-i + 2){l-z-i){l + z + d + i-2) 

{a{l + z + d)-2{2- 0) (2 + 2a{d + 1 + /) + a'^{d + 3)) 



{1 + z + d){l- z-^){l + z + d + ^-2) {I- z-2){l + z + d){l- z-^){l + z + d + ^-2) 
{a-imd + 3)i8-6^ + e) 



{l-z-2){l + z + d){l- z-Oil-z-^ + 2){l + z + d + ^-2) 

{a-imQ-5l + P){8-6^ + e) 

{l-z-2){l + z + d){l-z-Oil-z-^ + 2){l + z + d + ^-A){l + z + d + ^-2)' 

APPENDIX D: THE MATRIX P^K 

We defined the matrix K as 

^"'B'l^ir) = Jd'pe'^"^^, (Dl) 

where the elements are obtained by direct computation. Multiphcation with the transpose 
of the projector flB3p yields 

r (^-^) 

P^K = i2-' (D2) 

r(^^) ^ ^ 

where k is now a 2 x 4 matrix, 



z-l 



K = I ~ z+d+l {z+2-l)(z+d+l) z+d+l z+d+l 1 /-QgN 
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